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Abstract 

Exclusive electroproduction of <fi mesons has been studied in e ± p collisions at 
y/s = 318 GeV with the ZEUS detector at HERA using an integrated luminosity 
of 65.1 pb _1 . The 7*p cross section is presented in the kinematic range 2 < Q 2 < 
70 GeV 2 , 35 < W < 145 GeV and \t\ < 0.6 GeV 2 . The cross sections as functions 
of Q 2 , W, t and helicity angle 9h are compared to cross sections for other vector 
mesons. The ratios R of the cross sections for longitudinally and transversely 
polarized virtual photons are presented as functions of Q 2 and W. The data are 
also compared to predictions from theoretical models. 
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1 Introduction 



Exclusive electroproduction of vector mesons is a process which can be used to confront 
model predictions in a kinematic region that includes the transition between soft and hard 
dynamics. The process also lends itself to detailed experimental investigation; the decay 
products of the vector meson can be precisely measured so that distributions in several 
variables can be studied with high resolution over a large phase space. 

Many measurements of exclusive production of vector mesons, e p — > e V p, have been 
made at HERA [1-16]. In this paper, results from a study of exclusive <fi production are 
presented. The data sample represents a factor ~ 30 increase over previously published 
HERA results [9,10]. The measurements were made for virtuality of the exchanged photon, 
Q 2 , in the range 2 < Q 2 < 70 GeV 2 and for the photon-proton center-of-mass energy, W, 
in the range 35 < W < 145 GeV. The cross sections are presented as functions of Q 2 , W, 
the squared four-momentum transfer at the proton vertex, t, and helicity angle 0^. The 
W dependence was also extracted in bins of t, and the ratio, R, of the cross sections for 
longitudinally and transversely polarized virtual photons, determined from the angular 
distribution of the decay products of the (f> mesons, is presented as functions of W and 

Q 2 . 



2 Phenomenology 

A simple picture for the process e p — > e V p, where V represents a vector meson, can 
be formulated in the rest frame of the proton. In this frame, the virtual photon emitted 
from the electron fluctuates into a quark-antiquark dipole, which then scatters elastically 
off the proton. Long after the interaction, the qq pair forms a vector meson. It can be 
shown that, at high energy, the qq creation, scattering, and vector meson formation are 
well separated in time so that the cross section for the process can be factorized into 
terms representing the qq coupling to the photon, the dipole scattering cross section on 
the proton, and the final-state formation [17, 18]. The dipole scattering cross section 
on the proton depends on the transverse separation of the qq pair in the dipole. The 
interactions of large dipoles are thought to be primarily 'soft' and therefore described by 
Regge phenomenology [19]. On the other hand, the interactions of dipoles with small 
transverse separation of the qq pair are expected to be 'hard', and therefore calculable in 
perturbative QCD (pQCD). In lowest-order pQCD, the process proceeds via the exchange 
of two gluons between the proton and the qq dipole so that vector-meson production is 
sensitive to the gluon density in the proton [20]. 

The transverse size of the dipole is related to the transverse-energy scale of the interaction, 
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and can depend on Q as well as on the quark mass, M [21]. Large transverse-energy scales 
preferentially select small dipole sizes. The impact parameter for the dipole scattering 
on the proton depends on \f\t\. Large \t\ preferentially selects small impact parameters. 
The \t\ values in this analysis are small and no pQCD predictions are possible for the \t\ 
dependence in this regime. The pQCD predictions for the W and Q 2 dependence of the 
cross section should be more accurate as either or both of Q or M become large. 

The pQCD predictions depend on the square of the gluon density, leading to a steep 
dependence of the 7*p cross section on W, in contrast to expectations from Regge phe- 
nomenology. For example, a gluon density varying as xg(x) oc x -0 ' 2 , where x is the 
Bjorken scaling variable, would lead to <r(7*£> — > V p) oc W^ 0,8 . The Regge phenomenology 
expectation is <x(7*p — > V p) oc W ' 2 . The energy dependence of the 7*p cross section 
is therefore a good indicator of whether gluon exchange is dominant. Data from exclu- 
sive p production [2,4,5] show that the cross section cr(7*p — > pp) rises with W as W s , 
where 5 increases with Q 2 from about 0.2 at Q 2 = (photoproduction) to about 0.7 
at Q 2 « 30GeV 2 . However, in the case of exclusive J/ift production the cross section 
rises steeply with W even for photoproduction [12,13,16]. The investigation of elastic <fi 
production is of particular interest since the is an ss state, and the mass of the strange 
quark is not negligible in pQCD calculations. Experimentally, the extraction of the 
signal is very clean due to the narrow width of the state. 

The measurement of the variation of the W dependence of the cross section with t also 
provides a good test of the validity of pQCD calculations. In pQCD, little variation is 
expected, while for soft hadronic processes the W dependence varies as a oc W^ aF ®~ 1 \ 
with a P (t) = 1.08 + 0.25* [22,23]. 

Specific predictions for meson production are available in the dipole model realization 
(FS04) of Forshaw and Shaw [24]. The predictions depend on the strange-quark mass 
used in the photon wavefunction, on the parametrization of the dipole-proton scattering 
cross section, and on the vector-meson wavefunction. The gluon density does not appear 
explicitly but is implicit in the dipole scattering cross section. The predictions shown in 
this paper employed a strange-quark mass of 200 MeV, and a Gaussian wave-function for 
the 0-meson [25]. The dipole-proton cross section is obtained from the best fit to the 
total 7*p cross section [24] and implies a saturation of the cross section as x decreases. 

Specific predictions for the meson are also available in the model of Martin, Ryskin 
and Teubner (MRT) [26]. In the MRT model, the properties of the event are determined 
completely from the features of the photon wavefunction and the gluon density in the 
proton. Parton-hadron duality is invoked to relate the parton-level cross section to the 
cross section for the produced vector meson. 

As NLO corrections are not fully taken into account, the model calculations come with 
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significant normalization uncertainties. 

The data of this analysis are also compared to results from other vector mesons. The 
results are presented as a function of Q 2 + My, where My is the vector meson mass, for 
cross sections, and as a function of Q 2 /M v for the spin-density matrix element, r^, to 
test for scaling in these variables. 

A recent review of models and data can be consulted for more detailed information on 
vector meson production [27]. 

3 Experimental set-up 

The data were collected during 1998-2000 with the ZEUS detector and correspond to an 
integrated luminosity of 15.0 pb~ x for e~p and 50.1 pb _1 for e + p collisions with a proton 
energy of 920 GeV and an e ± energy of 27.5 GeV. Since no dependence on the lepton 
charge is expected, the two data sets were combined 1 . A detailed description of the ZEUS 
detector can be found elsewhere [28]. A brief outline of the components that are most 
relevant for this analysis is given below. 

Charged particles were reconstructed in the central tracking detector (CTD) [29] covering 
the polar-angle 2 region 15° < 6 < 164°. The transverse-momentum resolution for full- 
length tracks is cr(p T )/p T = 0.0058p T © 0.0065 © 0.0014/p T , with p T in GeV. 

The high-resolution uranium calorimeter (CAL) [30] consists of three parts: the forward 
(FCAL), the barrel (BCAL) and the rear (RCAL) calorimeters. Each part is subdivided 
transversely into towers and longitudinally into an electromagnetic section (EMC) and 
either one (RCAL) or two (FCAL and BCAL) hadronic sections. The CAL covers 99.7% of 
the total solid angle. Under test-beam conditions, the CAL single-particle relative energy 
resolution is a(E)/E = 0.18/v^E for electrons and a(E)/E = 0.35/^/E for hadrons, with 
E in GeV. 

The forward plug calorimeter (FPC) [31] was a lead-scintillator sandwich calorimeter with 
readout via wavelength-shifter fibers. It was installed in the beamhole of the FCAL and 
extended the pseudorapidity coverage of the forward calorimeter from rj < 4 to r\ < 5. It 
has since been removed to accomodate HERA magnets for the high-luminosity HERA II 
run. 



1 Hereafter, both e + and e are referred to as electrons, unless explicitly stated otherwise. 

2 The ZEUS coordinate system is a right-handed Cartesian system, with the Z axis pointing in the 
proton beam direction, referred to as the "forward direction" , and the X axis pointing left towards 
the center of HERA. The coordinate origin is at the nominal interaction point. 
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The small-angle rear tracking detector (SRTD) [32] consists of two planes of scintillator 
strips read out via optical fibers. It is attached to the front face of the RCAL and covers 
an angular range between 162° < 9 < 176°. The SRTD provides a transverse position 
resolution for the scattered electron of 0.3 cm [12] corresponding to an angular resolution 
of 2 mrad. 

The hadron-electron separator installed in the RCAL (RHES) consists of silicon diodes 
placed at a longitudinal depth of three radiation lengths. The RHES provides an electron 
position resolution of 0.5 cm if at least two adjacent pads are hit by the shower [33]. 

The luminosity was determined from the rate of the bremsstrahlung process ep — > e 7 p, 
where the photon was measured with a lead-scintillator calorimeter [34] placed in the 
HERA tunnel at Z = -107 m in the HERA tunnel. 



4 Kinematics and cross sections 

The following kinematic variables are used to describe exclusive <fi production, 

e(k)p(P) -> e(k') <f>(y) p(P') , 

where k, k', P, P' and v are, respectively, the four-momenta of the incident electron, 
scattered electron, incident proton, scattered proton and the 4> meson: 

• Q 2 = —q 2 = —{k — k') 2 , the negative four-momentum squared of the virtual photon; 

• W 2 = (q + P) 2 , the squared invariant mass of the photon-proton system; 

• y = (P ■ q)/(P ■ k), the fraction of the electron energy transferred to the proton in the 
proton rest frame; 

• x = Q 2 /(2P ■ q), the Bjorken variable; 

• t = (P — P') 2 , the squared four-momentum transfer at the proton vertex. 

The kinematic variables were reconstructed with the "constrained" method [2] which 
uses the momentum vector of the and the polar and azimuthal angles of the scattered 
electron. 

The ep cross section can be expressed in terms of the transverse, a^, and longitudinal, 
<7l, virtual photoproduction cross sections as 

-^T = r T (y,Q)(a T + ev L ), 

where IV is the flux of transverse virtual photons [35] and e is the ratio of longitudinal 
and transverse virtual-photon fluxes, given by e = 2(1 — y)/(l + (1 — y) 2 )- In the kinematic 
range studied here, e lies in the range 0.975 < e < 1, with an average value of 0.99. 
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The virtual photon-proton cross section, a 1 p ~*^ p = a T + ecr^, can be used to evaluate the 

'tot 



total exclusive cross section, (riot P = + &l, through the relation 



1 I p 

atot ~ TTTr 



where R = olI°t is the ratio of the cross sections for longitudinal and transverse photons. 
The helicity structure of (j) production is used to determine R as described in Section !^. 5. II 



5 Reconstruction and selection of events 

The signature of exclusive <p electroproduction, ep — > e <fi p, consists of the scattered 
electron and two oppositely charged kaons from the decay. The scattered proton is 
deflected through a small angle and escapes undetected down the beampipe. The data 
selection and analysis are described in detail elsewhere [36]. A brief outline is given here. 

The events were selected online by a three- level trigger [37,38]. The trigger required a 
scattered electron in the CAL with energy greater than 7 GeV, a minimum of two and a 
maximum of five tracks reconstructed with the CTD and less than 5 GeV in the towers 
of the FCAL closest to the beampipe. These cuts reduced the rate of background events 
while preserving high efficiency (> 99 %) for the signal events. 

The following criteria were applied offline to reconstruct and select the events: 

• the identification and energy measurement of the scattered electron used information 
from the CAL. The energy was required to satisfy E > 10 GeV. The impact point 
of the electron on the CAL was measured using three detectors: SRTD, HES and 
CAL. Given its superior position resolution, preference was given to the measurement 
from the SRTD. This was improved by the position obtained from HES or CAL when 
applicable. To ensure full containment of the electromagnetic shower and good position 
reconstruction, fiducial cuts were applied to the impact position of the electron on the 
face of the RCAL; 

• the mesons were reconstructed from the properties of the decay kaons. No particle 
identification was used. Events with two tracks of opposite charge each with pt > 
0.15 GeV and \r]\ < 1.7 were selected. These tracks were assigned the kaon mass 
and the invariant mass was formed. Track combinations with invariant masses falling 
within an allowed mass window were selected (see below). Events with additional 
tracks not associated with the scattered electron or with kaon decays were rejected; 

• the position of the reconstructed vertex was required to be compatible with that of an 
ep collision, |Zvtx| < 50 cm. The radial distance of the reconstructed vertex from the 
nominal beamline was required to be less than 0.8 cm to remove Kg — > 7i + tt~ decays; 
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• to remove events with large initial-state radiation, 45 < E — Pz < 65 GeV was im- 
posed, where the longitudinal energy-momentum variable E — Pz is calculated using 
the momenta of the two kaons and the scattered electron (the masses are neglected). 
This variable is peaked at twice the electron beam energy (55 GeV) for non-radiative 
DIS events. The E — Pz > 45 GeV cut removed events with a radiated photon of 
more than 5 GeV; 

• to suppress non-exclusive events, the energy of each CAL cluster not associated with 
either of the final-state kaons or the scattered electron was required to be less than 
0.3 GeV. To suppress further the contamination from proton-dissociative events, ep — ■> 
e <p Y, the energy in the FPC was required to be less than 1 GeV. These cuts restrict 
the mass of the proton-dissociative system, V, to My < 2.3 GeV. 

Events were required to be in a kinematic region where the properties of the final-state 
particles are well measured. Additionally, the kinematic range was limited to the region 
where the acceptance varies only slowly with the kinematic variables. This led to the 
following selection: 

Q 2 > 2 GeV 2 , 
\t\ < 0.6 GeV 2 , 
33.75 GeV + (1.25 GeV" 1 ■ Q 2 ) < W < 100 GeV + (3.7 GeV" 1 ■ Q 2 ) , 

with Q 2 given in GeV 2 in the last expression. The distribution of selected events in the 
x-Q 2 plane for 1.01 < m KK < 1.04 GeV is shown in Fig. [T] 

6 Monte Carlo simulation 

The acceptance and the effects of the detector response were determined using samples of 
Monte Carlo (MC) events. All generated events were passed through the standard ZEUS 
detector simulation, based on the Geant 3.13 program [39], the ZEUS trigger simulation 
package, and the ZEUS reconstruction software. 

The exclusive process ep — > e <p p was modelled using the Zeusvm [40] MC genera- 
tor interfaced to Heracles 4.6.1 [41] to account for first-order QED radiative effects. 
The parameters describing the W,Q 2 ,t and 6^ distributions in the MC simulation were 
adjusted such that the MC simulation reproduced the data distributions. The typical 
acceptance for -> K + K~ increases from 20% at Q 2 = 2 GeV 2 to 60% for Q 2 > 10 GeV 2 . 
The acceptance is small at the smaller Q 2 due to the fiducial cuts on the impact point of 
the electron. 
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Proton-dissociative events, ep — > e <fi Y , were modelled using the generator Epsoft [42]. 
The 7* p — * <p Y cross section was parametrized as 

with the same Q 2 and W dependence, f(Q 2 , W), as for the exclusive case. The values of 
b and f3 are discussed in Section 17721 



7 Extraction of the signal 

Figure El shows the invariant-mass distribution of the K + K~ pairs obtained after the 
selection described in Section A total of 3642 candidates were found after subtraction 
of non-resonant background in the mass range 1.01 < tukk < 1.04 GeV. The natural 
width of the resonance is comparable to the detector mass resolution, which is 1.8 MeV 
for tracks at central rapidity, deteriorating to ~ 5 MeV for tracks crossing the CTD 
through large angles. This resolution is well described by the MC simulation. 



7.1 Non-resonant background 

The non-resonant background was estimated for each bin in which a cross section was ex- 
tracted with an unbinned likelihood fit to the invariant-mass distribution. The assumed 
functional form was a p-wave relativistic Breit-Wigner convoluted with a Gaussian reso- 
lution function for the signal plus a background function with the shape a(rriKK — 2mx) b - 
The background, estimated from the fit, was subtracted from the number of observed 
events in the mass window. The contribution of the non-resonant background in the mass 
range of the signal is typically 18% at Q 2 = 2.4 GeV 2 , decreasing to 5% at Q 2 = 13 GeV 2 . 



7.2 Proton-dissociative background 

The remaining source of background consists of <p production accompanied by proton 
dissociation, e p — > e <fi Y, where the particles from the breakup of the proton are not 
detected. A similar study to that done for J ftp production [12] was performed. 

Proton-dissociative events with observed proton breakup were studied using a sample of 
diffractive events selected as described in Section with the following exceptions: 

• the elasticity criteria (last bullet in Section |3J) were not applied and the energy in FPC 
was required to satisfy -Efpc > 1 GeV; 
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• a pseudorapidity gap Arj > 3 was required between the energy deposits from the 
system Y and the decay products. 

The sample of proton-dissociation events identified with the FPC contained 300 <fi can- 
didates for \t\ < lGeV 2 in the kinematic range 45 < W < 135 GeV and Q 2 > 2GeV 2 . 
The W, Q 2 and 9^ dependences were found to be the same as for the elastic data sam- 
ple. The t dependence was found to have a slope (see Eq. 1) decreasing with Q 2 as 
b = (5.1 ± 1.3)/(1 + R(Q 2 )) GeV" 2 with R given in Section I8~5~T1 The MC distribution of 
My was tuned to describe the FPC energy distribution, yielding (3 = 2.0 ± 0.5. 

The fraction of proton-dissociative events in the final sample, averaged over t for \t\ < 
0.6 GeV 2 , was / p -di SS = (7.0 ± 0.4^2^)%; independent of W, Q 2 and 9h- The uncertainty 
is dominated by the uncertainty in the modelling of the My spectrum for My < 3 GeV, 
and by the simulation of the proton-remnant final state. The uncertainty of the behavior 
of the My spectrum at small masses was parametrized by allowing (3 to vary in the range 
1.5 < /3 < 3.0. The uncertainties on the simulation of the proton-remnant final state were 
estimated by a comparison of different MC simulations. 

The fraction of proton-dissociative events increases from 4.5% for < \t\ < 0.08 GeV 2 to 
14.5% for 0.35 < \t\ < 0.6 GeV 2 . The cross sections presented in the next sections were 
corrected for this background in bins of t, and globally in W, Q 2 and 9^. 



8 Results 



8.1 Cross sections 

In each bin of a kinematic variable, the ep cross section was extracted using the formula 

ep^e<f>p _ (j^data ~ j^bgdjCj- — fjp— diss) 



a 



ABC 



where -/Vd at a is the number of events in the data and iVb g d is the number of events from 
the non-resonant background. The overall acceptance is denoted as A, B accounts for the 
— > K + K~ decay branching ratio (49.2 ± 0.6)% [43], and C is the integrated luminosity. 

The total exclusive photon-proton cross section was calculated as 



\Q\W,t)=al^{Q\W,t). 



ep— »e <p p 
°DATA 

ep—>e cj) p 



In this way, the MC simulation was used to correct for radiative effects, detector accep- 
tance, cut efficiencies and the shape of the cross section within the bin. 



8 



The cross sections were measured for |i| < 0.6 GeV 2 . The results were then extrapolated 
to the full t range assuming da/dt oc e _b ''L The correction factor needed to extrapolate to 
the cross section integrated over t was evaluated bin-by-bin using the measured value of 
b for the given Q 2 , and ranged from 1.5 % to 5.0 %. The cross sections are quoted for the 
mass range 2m K < m KK < M + 5I\ where = 1019.4 MeV and T = 4.458 MeV [43]. 

8.2 Systematic uncertainties 

The systematic uncertainties on the measured cross sections were determined by varying 
the selection cuts and by modifying the analysis procedure. The sources of systematic 
uncertainties considered [36] were similar to those of previously published analyses [12,44]. 
Except for the two lowest Q 2 bins in the evaluation of da/dQ 2 , the systematic uncertainties 
(excluding normalization uncertainty) were smaller than the statistical errors. 

• The electron position reconstruction is critical in determining the acceptance and in 
the kinematic-variable reconstruction. Possible misalignments of the SRTD and CAL 
were estimated to be 2 mm, and the effect of such misalignments was tested via MC 
simulations. In addition, the electron fiducial cut was varied by 5 mm. The incomplete 
knowledge of the reconstructed electron position generally provided the largest source 
of uncertainty. 

• The elasticity cut is important in suppressing non-exclusive events. The energy dis- 
tribution of the most energetic cluster not assigned to the electron or one of the kaons 
was compared in data and MC simulation. The MC simulation was found to repro- 
duce the data distributions well, indicating a very small non-exclusive contribution 
with extra particles in the central region. The cut was varied from 0.3 to 0.4 GeV in 
both data and MC simulation to determine the uncertainties. 

• The stability of the fitting procedure and extraction of the signal was checked by 
reducing the fit range to rriKK < 1.07 GeV and changing the accepted mass range to 
1.01 < rriKK < 1-035 GeV, respectively. Except for bins with a small number of events 
the contribution to the overall uncertainty was small. 

• Uncertainties from the dependence on the MC parametrizations were also estimated 
and found to be small compared to other sources of systematic uncertainties. 

• The normalization uncertainty was different for different measurements. All mea- 
surements had a normalization uncertainty due to uncertainties in the integrated lu- 
minosity, ±2.5%, and due to the uncertainty in the decay branching ratio, ±1.2%, 
leading to a combined uncertainty of ±2.8%. In the case of the da/dQ 2 and da/dW 
measurements, the proton-dissociation background was subtracted globally and the 
normalization uncertainty increased to 1 5 ' 3 %. The uncertainty due to the subtraction 
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of proton-dissociative background for the da/dt measurement resulted in a t-dependent 
uncertainty as discussed in Section 17.21 For the angular distributions, the normaliza- 
tion uncertainty does not appear as only the shapes of the distributions were measured. 

The total systematic uncertainty (excluding the normalization uncertainty) was deter- 
mined by adding the individual contributions in quadrature. The correlated and un- 
correlated systematic uncertainties were evaluated separately. The typical size of the 
systematic uncertainty (excluding the normalization uncertainty) was 5% for the cross 
sections in W and Q 2 bins, and 5% at small \t\, increasing to 10% at \t\ = 0.6 GeV 2 . 



8.3 Dependence on W and Q 2 

The cross-section o~^ Q ^ ^, iiiccisiircd. els £i function of W and Q 2 , is shown in Fig. and 
given in Tabled I n addition to the measured cross section, this table gives the kinematic 
range over which the measurement was performed, the value of the kinematic variables at 
which the cross section is quoted, the acceptance, and the background-corrected number 
of <p events. The much higher statistics available for this analysis allow measurements of 
double-differential cross sections. 

The cross sections were fitted to a dependence a oc W s with results presented in Table 121 
The measured values of 5 show no Q 2 dependence within the present uncertainties. The 
values of S are compared to those from previous HI [4] and ZEUS [1,2,12,13] measurements 
in Fig. Eb, where the data are plotted as a function of Q 2 + M v . The values of 5 scale 
with this variable within the present uncertainties. 

The Q 2 dependence of the cross section, given in Table is shown in Fig. for W = 
75 GeV. The data are compared to previous ZEUS [8, 9] and HI [10] results. The <j) data 
from this analysis were fitted with a function of the form a oc (Q 2 + M?) _n expected 
in the VDM [45]. The new high-precision data show that the Q 2 dependence of the 
cross section cannot be fitted with a single value of n over the Q 2 range of this analysis. 
The fit parameter n was found to vary from n = 2.087 ± 0.055(stat.) ± 0.050(syst.) for 

2.4 < Q 2 < 9.2 GeV 2 to n = 2.75±0.13(stat.) ±0.07(syst.) for 9.2 < Q 2 < 70 GeV 2 . There 
are several possible causes that could lead to this behavior, including a variation of R with 
Q 2 , the dependence of as on Q 2 and the changing Q 2 dependence of the gluon density 
at fixed W. The longitudinal and transverse components of the cross section have been 
separately extracted using the measured value of R (see Section 18.5. 1|) . and are shown 
in Fig. EJd. The difference in the Q 2 dependence of the cross section for the two helicity 
components is clearly seen. 
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8.3.1 Comparison to Models 



The predictions of the MRT model [26] were compared to the data using three different 
gluon densities as shown in Fig. |Sk. The predicted W dependence with the ZEUS-S [46] 
and CTEQ6M [47] gluon densities are compatible with the data, whereas the predictions 
using the MRST99 [48] gluon density are too steep. The predicted Q 2 dependence is too 
steep for all gluon densities. 

The FS04 prediction [24] is also compared to the data in Fig. Efc,. The W dependence seen 
in the data is well reproduced, although the normalization is somewhat low at large Q 2 
and high at low Q 2 . A similar model [49] which did not employ a saturated dipole cross 
section showed a somewhat steeper W dependence. 

The expectations from the MRT and FS04 models for the Q 2 dependence are also com- 
pared to the data in Fig. ^jp. In this case, the ZEUS-S gluon density was chosen for the 
MRT model. This prediction agrees reasonably well with the data at higher Q 2 . At small 
Q 2 , the predicted transverse cross section is too high while the predicted longitudinal 
cross section is too low. The FS04 prediction reproduces the data better than the MRT 
prediction for the longitudinal component, and is very similar to the MRT prediction for 
the transverse component. 

8.4 Dependence on t 

The differential cross section, da2 ^^ p / d\t\, measured as a function of t in the range 
\t\ < 0.6 GeV 2 , is shown in Fig. for different values of Q 2 and W = 75 GeV. A function 
of the form da/dt = da/dt\t=o ■ e bt was fitted to the data and the results of the fit are 
given in Fig. Eb; along with measurements from other vector mesons [1,2,4,8,10,12,13]. 
The slope parameters from this analysis are given in Table EJ The values of b from 
this analysis show no Q 2 dependence within the present measurement uncertainties. The 
measurements in Fig. Et> are presented as a function of Q 2 + M v and are found to scale 
with this variable within the present data uncertainties. 

The data sample was analyzed to determine the W dependence as a function of t. In the 
Regge formalism, the differential cross section can be expressed as 

da/dt oc W A( - aF(f) - x \ (2) 

where ajp, the Pomeron trajectory, is usually parametrized as 

a P (t) = ajp(0) + a'jpt. (3) 

The trajectory measured in soft diffractive processes is ap = 1.08 + 0.25 t [22,23]. In 
contrast, a'jp is much smaller in J /if; photoproduction: ajp = 1.20 + 0.115 t [12]. In 
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this notation, ap = 1 + 5/4. The values of op(0) and Qp were determined by fitting 
the W dependence of the differential cross section at different \t\ values using Eq. (J2J) 
for Q 2 = 5 GeV 2 . Since the proton-dissociative process was found to have the same 
W dependence as the exclusive process, the extraction of a.p is not sensitive to this 
background contribution. The analysis was therefore extended up to |t| = 1 GeV 2 . The 
fit results are shown in Fig. UJi and are given in Table El The parameters of the trajectory 
were determined from a fit of Eq.© to the extracted ajp(t) values, as shown in Fig. [7)3. 

The results are: 

ajp(O) = 1.10 ± 0.02(stat.) ± 0.02(syst.); 
a'jp = 0.08 ± 0.09(stat.) ± 0.08(syst.) GeV" 2 . 

The value of a'jp is closer to the value measured in J / ip production than that measured 
in soft diffractive processes. 



8.5 Decay angular distributions 

The angular distributions of the decay of the <p provide information about the photon 
and </> polarization states. In the helicity frame [50], the production and decay of the <j) 
can be described in terms of three angles: &h, the angle between the <fi production plane 
and the lepton scattering plane; Oh, the polar angle, and <f>h, the azimuthal angle of the 
positively charged kaon. Under the assumption of s-channel helicity conservation (SCHC), 
the normalized angular distribution depends only on two angles, Oh. and iph — <Ph — ^h- 
The Oh distribution can be expressed in the form 

1 ^ ^l+r ° 4 + (l-3r ° 4 )cos 2 ^ . (4) 



N d cos Oh 



The spin-density matrix-element t[]q represents the probability that the </> is produced in 
the helicity-0 state from a virtual photon of helicity or 1. 

. ^ a i* 'p-*H> __ 

The normalized cross section - , tot la . at W = 90 GeV is shown in Fig.Blfor seven values 

a rfcos^y^) ° 1 

of Q 2 . The data were fitted to Eq. (jlj). The values of Tqq, determined from the fits, are 
given in Table |H1 and plotted as a function of Q 2 /My in Fig. El where they are compared 
to the values extracted for other vector mesons [1-3,8,9, 12, 13]. The values of r°o for the 
different vector mesons are found to scale in this variable. 

The values of Tqq were also extracted as a function of W for two different Q 2 values, and 
these are given in Table UJ No dependence on W was observed. 
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8.5.1 Longitudinal and transverse cross sections 



The ratio of the longitudinal to transverse cross section, R = (Tl/ct, wa s calculated as a 
function of Q 2 from according to the relation 

1 r 04 

R = l 7 M m> (5) 
e 1 - rgg 

which is valid under the assumption of SCHC 3 . The average value of e = 0.99 was used 
in extracting R. The values of R are presented as a function of Q 2 in Table El and as a 
function of W in Table [7| The scaling of Tqq with Q 2 /My implies that R also scales in 
this variable. This is expected in dipole models of diffraction [51]. 

The values of R are compared with previous ZEUS [8,9] and HI [10] results in Fig. ITUk as 
a function of Q 2 and in Fig.lTUb as a function of W. The Q 2 dependence is well described 
by the expression R = a(Q 2 / M^) b . The parameters were extracted from the fit to data 
for which the statistical errors are Gaussian, yielding a = 0.51 ± 0.07(stat.) ± 0.05(syst.) 
and b = 0.86 ± 0.11(stat.) ± 0.05(syst.). The prediction from the MRT model (using 
the ZEUS-S gluon density) is also shown in Fig. EH as is the FS04 prediction. The 
general power-law dependence of R with Q 2 is reproduced, but the model predictions 
systematically underestimate the measurements. The dipole-model prediction (FS04) is 
in somewhat better agreement with the ZEUS data than the MRT-model prediction (with 
the ZEUS-S gluon density), particularly at the lower Q 2 values. The weak dependence of 
R on W observed in Fig. ITUb is consistent with both the MRT and FS04 models. 



9 Summary 

The exclusive electroproduction of <fi mesons, e p — > e <fi p, has been measured with the 
ZEUS detector at HERA for photon virtualities in the range 2 < Q 2 < 70 GeV 2 , for 
photon-proton center-of-mass energies in the range 35 < W < 145 GeV and for four- 
momentum-transfer squared in the range |*| < 0.6 GeV 2 . The extracted 7* p cross section 
rises with W as a oc W s , with a slope parameter 5 ~ 0.4. This value is between the 'soft' 
diffraction value and the value observed in J/ip production. No Q 2 or t dependence of 5 
was observed within the present experimental precision of the data. 

The high-precision data from this analysis reveal that the Q 2 dependence of the cross 
section cannot be fitted with a single power over the measured Q 2 range. The longitu- 
dinal and transverse components of the cross section were separately extracted using the 

3 The validity of SCHC has been tested for other vector mesons and found to be adequate for the purpose 
of extracting R [3] . 



13 



measured value of R and the different Q 2 dependence of the cross section for the two 
helicity components is clearly seen. 

The t distribution, measured for \t\ < 0.6 GeV 2 , is well described by an exponential 
dependence over the range 2 < Q 2 < 70 GeV 2 . The slope ranges from 6.4 ± 0.4 GeV" 2 at 
Q 2 = 2.4 GeV 2 to 5.1 ± l.lGeV" 2 at Q 2 = 19.7 GeV 2 . 

The ratio of the cross sections for longitudinally and transversely polarized photons, R, 
increases with Q 2 and can be fitted by a power-law dependence. 

The data from this analysis were compared to previously published e p — > e V p data. 
The values of 5 and b were found to scale, within the present accuracy, when plotted as 
a function of Q 2 + M v . The ratio of longitudinal to transverse cross sections was seen to 
scale with Q 2 /M v . 

The MRT model does not reproduce the Q 2 dependence of the 7*p cross section observed in 
the data, while the FS04 prediction reasonably reproduces the data for the longitudinal 
photon polarization. This conclusion is also reflected in the better agreement of FS04 
with the data for the measurement of R vs Q 2 . The W dependence of the data can be 
reproduced in both models. 
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Q 2 range 


Q 2 


W range 


W 


Events 


A 


-y*p-*4>p 
w tot 


(GeV 2 ) 


(GeV 2 ) 


(GeV) 


(GeV) 




(%) 


(nb) 






35-45 


40 


203 


19.4 


76.4 ± 6.5±|;g 






45-55 


50 


255 


18.7 


101.2 ±7.7t^ 






55-65 


60 


220 


19.6 


101.9 ±8.3+^ 


2-3 


2.4 


65-75 


70 


210 


19.6 


112.8 ±9.4j£i 






75-85 


80 


167 


19.3 


107 illlg 






85-95 


90 


175 


20.7 


122 ± lltfo 






95-105 


100 


134 


19.5 


110 ± 111? 






40-50 


45 


184 


20.9 


47.0 ±4.1+^ 






50-60 


55 


146 


21.3 


44.3 ± 


3-5 


3.8 


60-70 


65 


169 


22.3 


56.7±5.1±t:o 






70-85 


77.5 


218 


22.0 


62.3 ± 5.0+^2 






85-100 


92.5 


158 


21.3 


57.4 ± 5.4+| 2 






100-115 


107.5 


123 


21.2 


59.0 ± 6.2t|° 






45-55 


50 


111 


34.6 


16.4 ± 1.8+1;} 






55-70 


62.5 


168 


37.6 


19.1 ± l.7t\l 


5-9 


6.5 


70-85 


77.5 


136 


38.5 


19.6 ± l-9±i;l 






85-100 


92.5 


123 


37.8 


21.6 ±2.3^;} 






100-115 


107.5 


116 


40.9 


23.li2.5tjj 






115-135 


125 


70 


44.6 


25.3 zt 3.5 ^ g 






50-60 


55 


64.5 


58.6 


5.05 ± 0.73+^ 






60-75 


67.5 


88.6 


57.0 


4.96 ± 0.59+°;^ 


9-30 


13.0 


75-90 


82.5 


87.3 


56.5 


6.12 ±0.81+°;^ 






90-105 


97.5 


86.6 


55.6 


6.82 ± 0.85lg;$ 






105-125 


115 


85.8 


58.0 


6.33 ± 0.78ig;i 






125-145 


135 


61.7 


56.8 


6.65 ± 0.971^ 



Table 1: The *y*p — > <\>p cross section, <rlJt~* > as a function ofW and Q 2 . The 
first uncertainty is statistical and the second systematic. The overall normalization 
error of t 5 ; 3 % is not included. The Q 2 and W values at which the cross section 
is evaluated are given in the second and fourth columns. The number of extracted 
events corrected for non-resonant background and the acceptance, A, are given for 
each measurement. 
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Q 2 range 
( GeV 2 ) 


Q 2 

( GeV 2 ) 


5 

(a oc W 5 ) 


2- 3 

3- 5 
5-9 
9-30 


2.4 
3.8 
6.5 
13.0 


0.41 ± O.lOlg;^ 
0.34 ± 0.13±g;g| 
0.43 ±0.15lg;^ 
0.38 ± 0.18±g;g^ 



Table 2: The S parameter (a oc W ) as a function of Q 2 . The first uncertainty 
is statistical and the second systematic. 



Q 2 range 


Q 2 


Events 


A 


u tot 


( GeV 2 ) 


(GeV 2 ) 




(%) 


(nb) 


2-3 


2.4 


1389 


19.3 


105.5 ± 3.4±|;g 


3-4.5 


3.6 


842 


21.4 


57.6 ± 2.4±|;§ 


4.5-6 


5.2 


420 


26.6 


31.1 ± l-8t\i 


6-8 


6.9 


376 


42.5 


17.9 ±i.i^;8 


8-11 


9.2 


314 


52.1 


11.06 ± 0.73±g;|f 


11-15 


12.6 


200 


59.9 


6.42 ± 0.52+°; 2 ^ 


15-20 


17.1 


61.8 


55.2 


2.50 ± 0.37lg;^ 


20-30 


24.0 


32.3 


62.5 


0.98 ±0.191^1 


30-70 


38.8 


10.5 


53.5 


0.37 ±0.1318$ 



Table 3: The *y*p — > <\>p cross section, 0~J o £~ , as a function of Q 2 for W = 
75 GeV . The second column gives the value of Q 2 at which the cross section is 
quoted. The number of extracted events corrected for non-resonant background 
and the acceptance, A, are given for each measurement. The first uncertainty is 
statistical and the second systematic. The overall normalization error o/l 5 ; 3 % is 
not included. 
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Q 2 range 


Q 2 


b ( oc e -b\t\) 

y a\t\ y 


da 
d\t\ 


t— n 


( GeV 2 ) 


( GeV 2 ) 


( GeV" 2 ) 


/~"u / ~\ 7"2\ 
(nb/ GeV J 


2-3 


2.4 


6 37 ± 32 + n nl 


678 ± 37 tf 6 


3-4.5 
4.5-6 

6-8 


3.6 
5.2 
6.9 


6.29 ± 0.42 toil 
5.26 ± 0.48 1°;^ 
5.49 ± 0.52 toiy 


352 ± 26 tH 
162 ± 16 t 9 13 
98 ± 10 t 8 10 


8-11 


9.2 


5.58 ± 0.54 toil 


61.3 ±6.6 tf 8 7 


11-15 
15-30 


12.6 
19.7 


5.45 ± 0.78 toll 
5.10 ± 0.98 t°oil 


32.9 ±5.5 tH 
8.8 ± 1.9 ti 8 



Table 4: T/ie s/ope parameter, b, as a function of Q 2 for W = 75 GeV . The cross 
section is quoted at the Q 2 value given in the second column. The first uncertainty 
is statistical and the second systematic. 



\t\ range 


|t| 


5 


a IP 


(GeV 2 ) 


(GeV 2 ) 


(a oc W s ) 


a IP = l + 5/ A 


0-0.08 


0.025 


0.34 ± 0.11±g;g| 


1.085 ± 0.0271^20 


0.08-0.20 


0.12 


0.39 ± 0.121°:°° 


1.098 ±0.030+°:° 22 


0.20-0.35 


0.25 


0.44 ±0.15^;H 


1.110 ±0.038l°;°28 


0.35-0.6 


0.45 


0.31 ±0.181^ 


1.078 ±0.045^086 


0.6-1.0 


0.73 


-on ±o.33i{j:t2 


0.972 ± 0.083^080 



Table 5: The 5 parameter and the Pomeron trajectory ajp as a function of \t\ 
The first uncertainty is statistical and the second systematic. 



Q 2 range 
( GeV 2 ) 



2-3 
3-4.5 
4.5-6 

6-8 
8-11 
11-15 
15-30 



Q 2 
GeV 2 ] 



2.4 

3.6 

5.2 

6.9 

9.2 

12.6 

19.7 



r 04, 

no 



0.529 ±0 
0.595 ±0 
0.680 ±0 
0.740 ± 
0.744 ± 
0.802 ±0 
0.825 ±0 



n9 cr+0.025 
uz,J -0.029 



031 



+0.036 
-0.030 



nqo+0.026 
• Uc >O-0.031 

n qo+0.028 
• UOO_ 027 

n4 o+0.017 
• u ^°-0.019 

n/1 o+0.022 
• u ^°-0.018 

Q fifi +0.019 
• uuu -0.024 



R = o L ja T 



1.13 
1.49 



0.12 
0.11 
+0.20 
0.18 



o i r+0.42 
z - id -0.33 

^•OO-0.50 

2-9418:" 

4.091J:S 



4.771?: 



+0.12 
-0.12 
+0.24 
-0.17 
+0.28 
-0.28 
+0.47 
-0.37 
+0.27 
-0.27 
+0.63 
-0.42 
49 +0.72 
58 -0.71 



The spin-density matrix- element rPi and the ratio of cross sections 



Table 6: 

for longitudinally and transversly polarized photons, R, as a function of Q 2 . The 
first uncertainty is statistical and the second systematic. Due to the transformation 
given in Eq. JSJ), the error on the measurement of R is asymmetric. 
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Q 2 range 


Q 2 


W range 


W 


r 04 
'00 




R = o L ja T 


(GeV 2 ) 


( GeV 2 ) 


(GeV) 


(GeV) 












35-50 


42.5 


0.597 ±0.038; 


H0.019 
-0.019 


1 cn+0.26 +0.12 
l.OU_ 22 _0.12 






50-65 


57.25 


0.569 ±0.038; 


(-0.017 
-0.017 


i OO+0.22 +0.10 
i -°°-0.19 -0.09 


2-5 


3 


65-80 


72.25 


0.511 ±0.041; 


hO.017 
-0.020 


i fiR+ - 19 + - 08 
1,uu ~0.16 -0.08 






80-95 


87.25 


0.611 ±0.043; 


1-0.026 
-0.028 


i cq+0.32 +0.19 
l.OJ_ 26 -0.18 






95-115 


102.75 


0.612 ±0.048; 


H0.033 
-0.033 


1 cq+0.37 +0.24 
l.OJ_ 29 -0.21 






45-60 


52.25 


0.800 ±0.042; 


1-0.031 
-0.029 


A nc;+l-33 +0.93 
^• UJ -0.87 -0.64 






60-80 


70.0 


0.767 ±0.038; 


hO.015 
-0.017 


o OO+0.84 +0.31 
°- oz -0.61 -0.30 


5-20 


8 


80-100 


90.0 


0.673 ±0.048; 


hO.015 
-0.014 


9 D8+ - 54 +0- 14 
z - uo -0.40 -0.13 






100-120 


110.0 


0.684 ± 0.050; 


1-0.032 
-0.030 


i q+0.59 +0.36 
^•J-y-0.43 -0.27 






120-145 


130.0 


0.839 ±0.074; 


H0.046 
-0.048 


e o+5.3 +2.5 
d -°-2.0 -1.4 



Table 7: T7ie spin- density matrix- element and i/ie raizo o/ cross sections for 
longitudinally and transversly polarized photons , R, as a function ofW for two bins 
in Q 2 . The first uncertainty is statistical and the second systematic. Due to the 
transformation given in Eq. ^Ej) the error on the measurement of R is asymmetric. 
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Figure 1: The distribution of the events in the kinematic plane of Bjorken x and 
Q 2 . The open grey dots represent all reconstructed events while the solid black dots 
are events in the accepted kinematic range (Section^. 
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ZEUS 




m KK (GeV) 

Figure 2: Invariant mass distribution of the K + K~ candidate pairs. The solid 
line shows the result of a fit including signal+background, while the dashed line 
shows the background component only. The shaded region indicates the mass range 
used for cross-section calculations. The error bars indicate the statistical uncer- 
tainties. 
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Figure 3: (a) Exclusive cross section as a function ofW for four values of Q 2 . 
The solid lines are the results of a fit to the form a oc W 5 . The overall normalization 
uncertainty of i 5 ' 3 % is not included in the error bar. (b ) The extracted values of 
5 compared with results from other vector mesons. The error bars represent the 
quadratic sum of the statistical and systematic uncertainties. Note that the ZEUS 
p data point at Q 2 + My near 8 GeV 2 has been shifted down by 0.3 GeV 2 for clarity 
of presentation. 
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Figure 4: Exclusive <fi cross section as a function of Q 2 + Mj for W = 75 GeV . 
(a) The total cross section compared with previous measurements, (b) The separate 
longitudinal and transverse contributions to the cross section. The curves represent 
the results of the fits to the total cross section as described in the text. The error 
bars represent the quadratic sum of the statistical and systematic uncertainties. The 
overall normalization uncertainty of + 5 ; 3 % is not shown. 
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Figure 5: (a) Exclusive <fi cross section as a function of W for three values 
of Q 2 compared with predictions from the MRT and FS04 model. Different gluon 
densities (ZEUS-S, MRST99, CTEQ6M) were employed in the MRT predictions 
as indicated in the figure. Note that the data at Q 2 = 3.8 GeV 2 are not shown 
in the figure for clarity, (b) Exclusive cross section as a function of Q 2 + M? 
for W = 75 GeV . The data are compared to the MRT and FS04 predictions. The 
ZEUS-S gluon density was used in the MRT model. 
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Figure 6: (a) Differential cross-sections do~/d\t\ for seven bins of Q 2 , for W = 
75 GeV and \t\ < 0.6 GeV 2 . The full lines are the results of fits to the form da/dt = 
do~/dt\t=o ■ e bt . The overall normalization uncertainty of 2.8 % is not included in 
the error bar. (b) The slope b, as a function of Q 2 + My, compared to other ZEUS 
and HI results. The inner error bars represent the statistical uncertainty; the outer 
bars are the statistical and systematic uncertainties added in quadrature. 
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Figure 7: (a) Differential cross-sections da/d\t\ as a function ofW for different 
values of t at Q 2 = 5 GeV 2 . The full lines are fits to the form W s . The mean 
contribution from proton dissociation has been subtracted. The overall normaliza- 
tion uncertainty of^\ % is not shown, (b) The extracted values of a p. The full 
line is a linear fit to ajp with the form given in Eq.ffl). The dotted line is the soft 
Pomeron trajectory [22]. The inner error bars represent the statistical uncertainty; 
the outer bars are the statistical and systematic uncertainties added in quadrature. 
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Figure 8: Normalized distributions of cos 9^ in seven Q 2 bins; the curves are the 
fits to Eq. (Op. The inner error bars represent the statistical uncertainty; the outer 
bars are the statistical and systematic uncertainties added in quadrature. 
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Figure 9: The values of r^ extracted from the fits shown in Fig. are plotted as 
a function of Q 2 /My and compared with values for other vector mesons. 
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Figure 10: (a) Ratio R = o~lIo~t as a function of Q 2 for exclusive <fi production; 
the full line is the result of the two-parameter fit as shown in the plot, (b ) Ratio 
R = OlI Ot as a function ofW extracted in this analysis for two different Q 2 values. 
The dashed curve in both plots is from the MRT model with the ZEUS-S gluon 
density, while the dotted curve is from the FS04 model. The inner error bars show 
the statistical uncertainty, the outer the statistical and systematic uncertainties 
added in quadrature. 
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